Critical Tsallis exponent in heavy ion reaction 
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The numerical solution of the nonlocal kinetic equation allows to simulate heavy 
ion reactions around Fermi energy. The expansion velocity and density profile show 
£NJ \ specific radial dependence which can be described with a Tsallis exponent of q = 5/3. 

This might be considered as an indication of a phase transition. 
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The question whether multifragmentation in heavy ion collisions around the 
Fermi energy is associated with a nuclear matter phase transition has been 
investigated for decades. From the Van der Waals equation of state it is ob- 
vious that, in order to describe such possible phase transition, we must have 



a kinetic equation capable to describe the second virial coefficient, or, more 
specifically, the excluded volume and the correlated pressure in equilibrium. 
^ ■ While the equilibrium virial correction including its quantum equivalent has 

been thoroughly investigated, it is astonishing that its nonequilibrium exten- 
sion has not been considered. Therefore, we have developed a nonlocal quan- 
tum kinetic theory which leads to an excluded volume as well as to a correlated 
pressure [1]. Solving the resulting kinetic equation demands no more numerical 
effort than solving the standard local Boltzmann equation extended by Pauli- 
blocking effects. Some experimental features of correlations are described by 
this result [2]. 

Searching for signals of a possible phase transition during the reaction is hin- 
dered by the over-shading of various reaction channels and four orders of 
magnitude later observation. Therefore dynamical models provide the correct 
approach to search for appropriate signals. In [3] it was found that in cen- 
tral heavy ion reactions around the Fermi energy a non-Hubblean expansion 
profile appears. This was associated with a peculiar density profile and it van- 
ishes towards standard Hubble expansion for higher energies. This has been 
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attributed to long range correlations which are typical for fluctuations near a 
phase transition. 

Since the heavy ion reaction is a complicated evolution of a correlated finite 
size system we might search for a simpler effective description of main fea- 
tures of the numerical solution. A promising short cut to describe effectively 
a statistics including finite size effects is the nonextensive statistical mechan- 
ics, starting from nonextensive entropies, as suggested by Druyvenstein [4], 
Renyi [5] , Sharma [6] or Tsallis [7] . For a discussion of the kinetics underlying 
generalized statistics see [8]. Here it will be shown that both the velocity as 
well as the density profile can be associated with an anomalous diffusion of 
fractional derivative Fokker- Planck- like equation [9] with a Tsallis exponent 
q = 5/3. Since the latter value represents the border between GauB and Levy- 
like fluctuations [10], we consider it a hint of phase transition. 

Within the (extended) quasi-particle and quasi-classical approximations [11] 
we keep the gradient terms in the scattering integral of the Kadanoff and 
Baym equation and obtain a quantum non-local kinetic equation [11], 



The arguments of the distribution functions and the corresponding ones of the 
quasiparticle energies e are nonlocal, /i = f(k,r,t), / 2 ~ = /G°? r— A 2 , t), f% = 
f(k-q-A K ,r-A 3 ,t-A t ), and f 4 = f(p+q-A K ,r-A 4 ,t-A t ). The scattering 
probability is the square of the T-matrix, P~ = j0jsb (ei+ej — — ej — 2A.e) 

T (e 1+e --A E ,k-^,p-^,q,r-A r ,t-f) 
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. All non-local corrections 



are given by derivatives of the scattering phase shift = Im lnT(f2, k, p, q, r, t) 



. d<p 1 d(f> 18(f) 

2 dp dq dk' 3 dk ) 4 dk dq 

The collision is of finite duration A t . During this time particles can gain mo- 
mentum and energy Ak,e due to the medium effect on the collision. Three 
displacements A2,3,4 correspond to the initial and final positions of two collid- 
ing particles/holes. The numerical values of the shifts calculated with realistic 
potentials are available [12] and lead to simple off-set in the the algorithm 
which simulates the collision [2]. 

We plot in Fig. 1 the time dependence of the power of radial velocity v oc 
r a for different bombarding energies and the corresponding power of radial 
density dependence in Fig. 1. There are two distinct behaviors dependent 
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Fig. 1. The time dependence of the power law velocity profile (left) with respect to 
the radius v oc r a for different lab energies [3]. The dotted lines show the surface 
matter behavior and the solid lines depict the bulk matter behavior which was 
distinguished at the radius R ~ lOfm. The time dependence of the power law fit of 
the particle density n oc r~@ in the surface region (R > 10 fm) is seen on the right 
hand side plot. 



on the bombarding energy. For higher energies, 60MeV and 90MeV, we see 
that the velocity exponent approaches a ~ 1, in agreement with the Hubble 
expansion. This is associated with a vanishing density exponent. The slightly 
negative values in Fig. 1 can be understood as ring-like expansions of matter 
(pancakes). The interesting observation is that for energies below or around 
Fermi energy the density exponent is (3 ~ 3 and the surface velocity exponent 
a ~ 2. This remarkable feature is present in very different simulations we 
made with different reactions. It is more clearly pronounced in the nonlocal 
scenario when compared to the local scenario. Therefore, we propose that it 
is connected with correlations. 



Both exponents can be understood from a critical Tsallis exponent q = 5/3. 
To see this, we first interpret the radial and time dependent density as a prob- 
ability distribution in the sense of a Fokker-Planck equation with fractional 
derivatives, the solution of which reads [9] 
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with z = x(\k\\t) t 2 -t+ x . The order of fractional derivative 7 is linked to the 
Tsallis exponent by q = ^| ~ |. For a critical value of 7 = 2 or q = 5/3, one 
obtains indeed the above observed cubic density profile P 2 {x 1 t 1 ^) oc x~ 3 . The 
associated velocity profile can be explained also by observing that for 7 = 2 
from (3) follows 
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and the velocity profile can be estimated to 

1 1 1 z 2 

v(r) ~< xP 2 (x, t) >oc ^(H 1 + bz) + — ) « — + — + o(bz 3 ). (5) 

tr 1 + bz b z 2 

This explains the radial quadratic exponent of the observed velocity profile. 

In conclusion, we saw that for bombarding energies below the Fermi energy 
there appears a peculiar quadratic velocity profile with respect to the radius 
and an associated cubic power law of the density Both features can be fit- 
ted by a fractional derivative Fokker-Planck equation with the critical Tsallis 
exponent q = 5/3. Since this exponent marks exactly the transition between 
Gaufi- and Levy-like diffusion we suggest that the observed radial dependence 
might be associated with a phase transition during the reaction stage. 
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